What Determines Inhomogeneous Line Widths in Semiconductor Microcavities? 
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Microcavity inhomogeneous line widths are calculated numerically using a microscopic two dimensional model of the 
polariton interaction with quantum well disorder. The calculations show show that in most structures the line widths 
are determined by disorder scattering between polariton and exciton states. This is because motional narrowing 
effectively removes the contribution due to multiple scattering between polariton states. 



Recent experimental results ^ have shown the im- 
portance of motional narrowing in determining the line 
widths of features in the reflectivity spectra of semicon- 
ductor microcavities. In a microcavity, the fundamental 
optical excitations are polaritons, arising from the cou- 
pling of a confined cavity photon mode to exciton states 
in quantum wells within the cavity. The spectral fea- 
tures associated with the polariton are broadened, both 
homogeneously due to the finite life time of the cavity 
photon, and inhomogeneously due to the interaction with 
disorder in the quantum wells. Motional narrowing oc- 
curs because the polaritons, being part photon, have very 
long wavelengths compared with those typical for exci- 
tons. The inhomogeneous broadening is therefore sub- 
stantially reduced by averaging over a large area of the 
much shorter length scale disorder potential. 

The important physical parameters which determine 
the inhomogeneous line widths arc the dispersions of the 
exciton and photon, the strength of their interaction, and 
the statistical properties of the disorder. A good ap- 
proximation is to treat each dispersion as parabolic, the 
photon mass, ~ 3 x 10~^me, being a consequence 
of the effects of cavity confinement The exciton- 

photon coupling leads to the formation of two polariton 
branches, separated at resonance by the vacuum Rabi 
splitting, hfl ~ 5 meV, which provides a measure of the 
strength of the interaction. Finally, the quantum well 
disorder potential is characterised by its amplitude Vq, 
of order a few meV, and a correlation length 1^ 100 A. 

The original treatment of motional narrowing in Rcf. 
assumed the weak disorder limit, Vo ^ hQ, in which 
only scattering between low momentum polariton states 
on the same branch is allowed, ft was then possible to 
use a simple scaling argument to predict the variation 
of the line widths with detuning. However, in a com- 
ment on that model, Agranovich et al. showed that 
the same scaling argument predicts the actual value of 
the inhomogeneous line width, at resonance, to be 



r2d ^ 



\2Mill 



(1) 



For physically reasonable parameters, is extremely 
small, ^ 10^^ meV, compared with the experimental in- 
homogeneous line widths of about 0.3 meV. 

The main aim of this letter is to demonstrate that, 
besides the polariton multiple scattering processes which 
give Eq.(Q), there is a contribution to the line width due 



to disorder scattering between polariton and higher mo- 
mentum exciton states. Such scattering is not included 
in the scaling treatment, since it assumes the polariton 
dispersion is parabolic, while in fact the two branches 
decouple at large wave-vectors. The more realistic cal- 
culation described here uses a numerical simulation of a 
polariton in the presence of the disorder, which automati- 
cally includes both types of scattering. However, in a typ- 
ical structure, the small value of the polariton multiple 
scattering contribution, estimated above, means that po- 
lariton to exciton scattering processes dominate. On the 
other hand, in a high quality sample, where Vq < 
scattering of polaritons into exciton states becomes en- 
ergetically impossible, and only the small polariton mul- 
tiple scattering contribution remains. 

The numerical model used here is physically similar, 
though solved differently, to the one dimensional simula- 
tions of Savona et al. However, motional narrowing 
in one dimension is significantly less important than in 
two, because averaging over a Id strip samples much less 
of the variation of the potential than an average over a 2d 
area of equivalent size. This difference can be seen very 
clearly by calculating the Id equivalent of Eq.(^, for the 
contribution of polariton multiple scattering to the line 
width. 
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For the parameters that gave ~ lO^"' meV, Fid ~ 
10~^ meV, three orders of magnitude greater, and far 
more comparable to the experimental line widths. This 
comparison suggests that the balance between polariton 
multiple scattering and polariton to exciton scattering 
should be very different in such Id calculations when 
compared to real 2d structures. 

The numerical treatment is based on the two level 
model of the polariton, in which a discrete cavity mode 
couples to the exciton ground state. The exciton inter- 
acts with the quantum well disorder potential, so in-plane 
momentum is not conserved. Thus the model consists 
of coupled two dimensional partial differential equations 
describing the in-plane motion of both the exciton and 
cavity photon. This model contains a number of approx- 
imations, most notably the omission of excited exciton 
states and the lack of a proper treatment of the effects 
of the magnetic field applied in the experiments. These 
approximations restrict the possibilities for precise com- 
parisons with experiment, but the the good general agree- 
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ment which is obtained suggests that the essential physics 
is included. 

The Hamiltonian for the model system takes the form 
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In addition to the photon mass Mi , and the vacuum Rabi 
splitting, Ml^ discussed above, Mg is the exciton mass, 
5 the detuning of the cavity photon relative to the exci- 
ton, and 7, the photon homogeneous width (due to es- 
cape through the mirrors). The quantum well disorder, 
V{y), is a Gaussian stochastic potential, constructed in 
the Fourier domain, following the method described by 
Glutsch et al. ||]. Its amplitude, Vq, is defined as the 
width at half maximum of the potential probability dis- 
tribution. 

The inhomogeneously broadened spectra are calcu- 
lated from the photon Green's function 

Gfc(i) = -i< tiu\e-'""''\tik > e{t) (4) 

where > represents a state of the system (not an 
eigenstate) consisting of a plane wave cavity photon with 
in-plane wavevector k and no exciton. All the results 
given in this letter are for normal incidence, correspond- 
ing to fc = 0. The Green's function describes how a pho- 
ton which enters the cavity is scattered out of its initial 
plane wave state by interactions with the disordered exci- 
ton. The plotted spectra are actually the photon spectral 
functions, obtained by Fourier transforming the Green's 
function to give Gfe(tj) and extracting the imaginary part. 
They correspond to what would be seen in an absorption 
measurement, but the line widths should be the same in 
reflectivity experiments. 

The Green's function is obtained using an approach 
similar to that of Glutsch et al. Starting from an 
initial state |V'(0) >= \f^k >, the wavefunction \ip{t) > 
is calculated by numerically solving the time depen- 
dent Schrodinger equation using the Hamiltonian Eq.(^. 
From the solution at each time t, Gk{t) is found by evalu- 
ating < ^fclV'(i) > = < Mfc|e~*^*^''|Aife >■ The calculation 
is carried out on a two dimensional spatial grid, using 
a standard alternating direction Crank-Nicholson algo- 
rithm 1^. The main difficulty is the need to use a grid 
which is fine enough to reveal the correlated structure 
of the potential, on a length scale Ic ~ 100 A, yet large 
enough to avoid significant size quantisation for the low 
mass photon states. The present numerical results are 
for a 2^^ X 2^^ grid with spacing 100 A, giving photon 
quantisation energies 0.05 meV, considerably smaller 
than typical calculated line widths of ^ 0.3 meV. 

Fig.|l| shows spectra calculated for a realistic set of pa- 
rameters: M; = 3 X 10~^me, Me = 0.5me, Tifi = 5 meV, 
Vo = 3.5 meV and Ic = lOOA. For clarity, the photon ho- 
mogeneous line width, 7, has been set to zero. The spec- 
tra display the variations in intensity typical of an anti- 
crossing between two states, only one of which, the pho- 
ton, has intrinsic strength. Close to resonance, where the 
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FIG. 1. Numerical absorption spectra (dots) at various val- 
ues of the detuning 5. The dashed curve shows the bare ex- 
citon lineshape. The solid lines are the predictions of the 
coupled oscillator model (see text). 

polariton branches are equal mixtures of photon and exci- 
ton, the two peaks have similar integrated strength. Fur- 
ther away from resonance, the more photon-like peak is 
strong, while the other, predominantly exciton, is weak. 
The figure also shows the bare exciton spectrum, which 
is asymmetrically broadened because of the finite exciton 
mass f7|, |l. 

In Fig^g, the line widths measured from the spectra 
are plotted as a function of detuning. The asymmetry 
between the line widths of the two branches at zero de- 
tuning is very apparent: the upper branch has a width 
of ~ 1 meV, and the lower branch 0.25 meV. When a re- 
alistic photon homogeneous broadening of 7 = 1.25 meV 
is included in the calculation, the widths become 1.4 and 
0.6 meV respectively, close to the experimental values 
of ^ 1.5 and 0.75 meV found in Ref. [|l|. However, be- 
cause of the lack of a proper treatment of the magnetic 
field used in the experiments, this agreement can only be 
regarded as suggestive. 

The solid lines in Figs.|l|j^ show the results of a simpler 
model, which provides an excellent fit to the numerical 
data. In this model, the full exciton Green's function 
G'f.f.,{uj) is approximated by its diagonal part 5kk'G'j.{uj). 
With k thus conserved, the polariton problem simplifies 
to a pair of coupled oscillators, which is easily solved to 
obtain the polariton Green's function 
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where Sk = 5 + ?i^fc^/(2M/) is the detuning for wave- 
vector k. The bare exciton Green's functions, G^(a;) are 
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FIG. 2. Numerical line widths as a function of detuning, 
5, for upper (open symbols) and lower (filled symbols) polari- 
tons. The solid lines are the predictions of the coupled oscil- 
lator model, the dashed lines those of the absorption model. 

calculated using the same type of numerical simulation as 
for the full polariton. However approximate expressions 
for the finite mass exciton spectral function, by Glutsch 
and Bechstedt could be used for large exciton masses. 

Eq. excludes polariton multiple scattering, since the 
polariton always has a definite wave-vector, k. However, 
polariton to exciton scattering processes are included, be- 
cause the exciton Green's function G%{i^) takes care of 
the multiple scattering of cxcitons between different mo- 
mentum states. It has a finite imaginary self energy de- 
termining the rate at which excitons scatter out of the 
state with wave- vector k. The good fit to the numerical 
data which is obtained indicates that the contribution of 
polariton multiple scattering is negligible for these pa- 
rameters, as the estimate of Eq.(|l|) suggests. 

In the coupled oscillator model, the difference between 
the widths of the two polariton branches is a consequence 
of the asymmetry of the exciton line shape, which is, in 
turn, a result of the finite exciton mass. On the low 
energy side of the line, there is an exponential cut-off, 
reflecting the distribution of minima in the disorder po- 
tential, while on the high energy side the strength falls off 
more slowly, as ^ determined by the perturbation 
of higher momentum plane wave states by the potential 
. This explanation for the asymmetry is thus similar 
in essence to the suggestion by Savona et al ^ that the 
larger upper branch width is caused by scattering into 
higher momentum exciton states. 

An approximate expression for the polariton line 
widths can be found by examining Eq. (^) . For sufficiently 
small broadening, the polariton lineshape approximates 
to a Lorentzian, with full width F dependent on the imag- 
inary part of the denominator at the polariton energy, ujp, 
according to 



F = 2|Qp(7-f}V4Ini{G|(^p)}) (6) 

where |c/p is the photon fraction of the polariton. This 
approximation is shown as the dashed lines on Figj^. 
For large photon fractions, when the polaritons are res- 
onant with the tails of the exciton line (see Fig.|^), the 
approximate expression agrees fairly well with the ex- 
act curve. However, when |qP becomes small, Eq.(H) 
predicts F ^ 0, while the true value tends to the exci- 
ton line width. This behaviour can be understood more 
physically by noting that the second term in Eq.(|^) de- 
scribes homogeneous broadening of the polariton due to 
the loss of photons in the cavity by absorption into ex- 
citon states - it is proportional to the photon fraction 
and the exciton absorption strength — Im{G^.(wp)}. This 
is to be be expected for a state in the tail of the ex- 
citon line shape, which is too weak to require a strong 
coupling treatment, and thus acts simply as a source of 
absorption, like the continuum states without a magnetic 
field 1^. The absorption picture breaks down when the 
polariton is resonant with the states near the centre of 
the exciton line, and the full strong coupling treatment 
is required. 
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FIG. 3. Numerical line widths as a function of disorder 
strength, Vb, for upper (open symbols) and lower (filled sym- 
bols) polaritons. The solid lines are the predictions of the 
coupled oscillator model. 

The coupled oscillator model described above can be 
shown to be a generalisation of the phenomenological 
treatment of a translationally invariant Gaussian broad- 
ened exciton, proposed by Houdre et al |jll[. The main 
improvement is that the Gaussian exciton line is replaced 
by an asymmetic form because of the finite exciton mass. 
More importantly, the present work demonstrates that 
motional narrowing is the reason why such a model suc- 
ceeds, in that it eliminates the polariton multiple scat- 
tering contribution to the line width. 

The most intriguing feature of the Gaussian model, 
discussed by Savona and Weisbuch , is the prediction 
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that the disorder contribution to the hne width should 
disappear in high quahty structures when polariton to ex- 
citon scattering becomes impossible. In terms of the ab- 
sorption picture, this happens when the polariton states 
are so far into the tails of the exciton line that no absorp- 
tion occurs. The present calculations support this predic- 
tion, as is shown by Fig.^, where numerical and coupled 
oscillator model line widths are plotted as a function of 
disorder strength, Vo- As in Rcf. jl^, the lower branch 
width rapidly becomes very small when Vq < n/2. The 
effect is much less pronounced for the upper branch, be- 
cause of the longer tail on the high energy side of the 
finite mass exciton line shape. 
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FIG. 4. Numerical line widths as a function of photon 
mass, Ml, for upper (open symbols) and lower (filled sym- 
bols) polaritons in (a) two dimensions and (b) one dimension, 
with Me — > oo. Also shown (dashed lines and right hand 
scale) are the polariton multiple scattering expressions, 
and Fid of Eqs. (|l|) and (1). 

All the numerical data in Figsj^^ are very well ex- 
plained by the coupled oscillator model, which includes 
only polariton to exciton scattering. It is not possi- 
ble in the numerical calculations to resolve directly the 
very small polariton multiple scattering term, from 
Eq.(|l]), even for Vq ^ fi/2, when it should be the only 
contribution to the line width. However, indirect evi- 
dence for the existence of the process can be obtained 
by considering non-physical situations in which polariton 
multiple scattering is more important. As an example, 
Fig.^(a) shows the dependence on the photon mass. Mi, 
of the numerically simulated line widths at resonance. 
The oscillator model predicts that the line widths of the 
two branches should be equal, as the calculation is for 
infinite exciton mass, and independent of Mi. For the 
real Mi = 3 x lO^^me, the line widths are indeed equal. 



However, when Mi is significantly increased a difference 
appears, which is clearly due to polariton multiple scat- 
tering, as it grows in proportion to F2d, plotted as the 
dashed curve on the figure. Fig.^(b) shows similar data 
for a one dimensional system. Again, the difference be- 
tween the upper and lower branch widths grows, for small 
Ml, in proportion to Fid ||l3|| . However, because motional 
narrowing is much less effective in Id than 2d, even for 
the real M; = 3 x 10~^TOe, there is a measurable dif- 
ference in line widths, some 0.05 meV. This result em- 
phasises that there are significant differences between the 
scattering processes which occur in Id and 2d systems. 

In summary, the numerical simulations presented in 
this letter have shown that there are two types of pro- 
cesses by which disorder contributes to the broadening of 
microcavity line widths. In typical structures, scattering 
between photon and exciton states dominates. From the 
point of view of the photon, this is just a form of life time 
broadening, caused by absorption into states in the tails 
of the exciton distribution, and is thus in some senses ho- 
mogeneous. The true inhomogeneous width, the result of 
multiple scattering of polaritons by the disorder, is ex- 
tremely small because of motional narrowing. However, 
it should be the only disorder contribution to the width 
of the lower polariton branch in high quality structures. 
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this fits better to calculated exciton line widths in Id. 



4 



